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TEST PARTICLE PROPAGATION IN 


MAGNETOSTATIC TURBULENCE 
m. THE APPROACH TO EQUILIBRIUM 

A, J. Klimas 
G. Sandri 
J. D. Scudder 
D. R. Howell 

ABSTRACT 

The asymptotic behavior, for large time, of the quasi-linear diabatic 
solutions and their local approximations is considered. A time averaging pro- 
cedure is introduced which yields the averages of these solutions over time 
intervals which contain only large time values, A discussion of the quasi-linear 
diabatic solutions which is limited to those solutions that are bounded from be- 
low as functions of time is given. It is shown that as the upper limit of the time 
averaging inte.-’val is allowed to approach infinity the time averaged quasi-linear 
diabatic solutions must approach isotropy (,ii -independence). The first deriva- 
tive with respect to ,u of these solutions is also considered. This discussion is 
limited to first derivatives which are bounded functions of time. It is shown 
that as the upper limit of the time averaging interval is allowed to approach 
infinity the time averaged first derivative must approach zero everywhere in jj. 
except at ^ = 0 where it must approach a large value which is calculated. The 
impact of this large derivative on the quasi-linear expansion scheme is dis- 
cussed. An H- theorem for the first local approximation to the quasi-linear 
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diabatic solutions is constructed. Without time averaging, the H-theorem is 
used to determine sufficient conditions for the first local approximate solutions 
to asymptote, with increasing time, to exactly the same final state which the 
time averaged quasi-linear diabatic solutions must approach as discussed 
above. An explanation, based on our numerical results, of how this peculiar 
final state develops is given. A formal discussion of the integro-differentiai 
equation, equation n.l of pap .r II. in this series, which contains all of the funda- 
mental approaches to statistical mechanics is given. The exact solution of this 
equation, as well as all of its local approximations is considered. For this dis- 
cussion the local approximation concept is generalized to include all of the 
approximations which result when any of the Fade' approximants to the Laplace 
transformed kernal of equation n.l are substituted for the exact kernel trans- 
form instead of only those Fade' approximants which lie near the diagonal of the 
Fade' table and correspond to the kernel convergents. This discussion is limited 
to those solutions and their local approximations which ai*e bounded from below 
as functions of time, and to those local approximations which are generated by 
a Fade' approximant whose range is not reduced compared to the range of the 
exact kernel transform. Then, it is shown that the time average of the exact 
solution and also of all its local approximations must asymptote, with increasing 


time, to a unique final state which is calculated. 




TEST PARTICLE PROPAGATION IN MAGNETOSTATIC TURBULENCE 


III. THE APPROACH TO EQUILIBRIUM 
I. INTRODUCTION 

In the first paper of this series^ (paper I.} we demonstrated the failure of 
the familiar Fokker- Planck description of the propagation of a charged test 
particle in a simple model of plasma turbulence. It had been shown previously^ 
that in the axi- symmetric slab model of the magnetostatic turbulence the 
Fokker-Planck equation for tlie ensemble and gyro-phase averaged probability 
distribution fimction follows from the quasi-linear diabatic equation through the 
application of the adiabatic approximation. We showed in paper I. that the adia- 
batic approximation is the leading term in the standard linear time-scale^ ex- 
pansion in a small parameter. As is typical in the derivation of many kinetic 
theories this small parameter can be expressed as the ratio of two time scales 
which are relevant to the mechanical system under consideration and which are 
assumed to be well separated; i.e., one of the time scales is assumed to be much 
larger than the other so that an appropriate ratio can be formed which is small 
and which then plays the role of an expansion parameter. In paper I. we showed 
that for the test particle propagation problem the assumed separation of the two 
relevant time scales is generally satisfied if the plasma turbulence is weal^, 
however, if a strong mean magnetic field is present there is a small but impor- 
tant region of the particle phase space in which the two time scales become 
comparable, or mixed. This mixing of the two time scales was shown to lead to 
a non-uniformity in the expansion which contains the Fokker-Planck result as 
the leading term. This failure occurs in the vicinity of 90“ pitch angle and is 
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due to the presence of a strong mean magnetic field. Thus, we showed that the 
failure of the Fokker- Planck description is due to the mixing of time scales 
which were incorrectly assumed well separated. 

In the second paper in this series'* (paper II.) we introduced our new local 
approximation method. In that paper we showed that this approximation method 
can be applied to the fundamental equations of statistical mechanics if they are 
cast in the integrodifferential form of equation II. 1 . We demonstrated the use of 
this new method through its application to a very simple exactly soluble model 
equation and to the quasi-linear diabatic equation whose solution, as discussed 
above, could not be well approximated through the use of the more standard 
approximation techniques. We showed that for these problems the local approxi- 
mation method does not require a small parameter for expansion purposes. 
Instead, the equation under consideration was Laplace transformed and approxi- 
mated near the origin of the Laplace variable using the smallness of the Laplace 
variable in that vicinity. The purpose of this particular approach is identical to 
that of the more standard approximation methods; i.e., we seek an accurate 
representation of the slow and/or long time evolution of the probability distribu- 
tion function. However, in the local approximation method, since we do not re- 
quire a small parameter, we avoid the non-uniformity problems associated with 
the standard asymptotic expansion procedures. In particular, this method was 
used to replace the imsuccessful Foldcer- Planck approximation with the re- 
markably successful first local approximation. 

Except for our considerations of the very simple exactly soluble model 
equation in paper I., our characterizations of success or failure of the approxi- 
mation methods discussed above rest on the results of numerical integrations 
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of the equations which generate the quasi-linear diabatic solutions, the quasi- 
linear adiabatic (Fokker- Planck) solutions, and the first local solutions. The 
lengths of these numerical integrations are necessarily limited, and therefore, 
our comparisons of the approximate solutions with the exact ones typically ex- 
tend over two to three Larmor periods of elapsed time. While this amount of 
time is sufficient to reject the Fokker- Planck approximation we are still left 
with many unanswered questions concerning both the quasi-linear diabatic solu- 
tions and their first local approximations. Do the diabatic solutions approach a 
final steady state? If so, what is this state? Do the first local approximations 
to the diabatic solutions approach a final steady state also? If so, how does this 
final state compare with that of the diabatic solutions? These questions are 
answered in this paper. 

We iind that the accuracy of the first local approximations to the diabatic 
solutions, which is evident during the first few Larmor periods from the 
numerical results, almost certainly persists for ^ time beyond the origin 
where the initial conditions on these solutions are specified. 

In section tj, of this paper we discuss the behavior of the quasi-linear 
diabatic solutions for very large times. We limit our discussion to solutions 
which are bounded from below as functions of time and consider magnetostatic 
turbulence whose power spectrum contains finite total power and has no zeroes 
in its high frequency tail except at infinite frequency. With these limitations we 
are able to show that the time average of these solutions, as defined by equa- 
tion 111,1 below, must asymptote, with increasing time, to a final state which is 
isotropic (independent of >i). Even so, we are also able to show that the first 
derivative with respect to q of these solutions, similarly time averaged, must 
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asymptote to a large non-zero value at aj. = 0. In the last part of this section 
we malce the further assumption that the solutions under consideration have a 
first derivative with respect to ^j. which is a bounded function of time. Then we 
are able to show that this derivative, again time averaged, must asymptote to 
zero everywhere that jx ^ 0. Various evidence from our numerical results is 
presented later in this paper to give an indication of how this peculiar final state 
actually develops. 

We have noted earlier (in paper I.) that the structure of the quasi-linear 
diabatic equation does not allow us to exclude the possibility that its solutions 
contain wave-like behavior for large times. Thus, the solutions to the quasi- 
linear diabatic equation need not approach a final state even though their time 
averages must. On the other hand, we have never found evidence from our 
numerical solutions, for as far as they go in time, for any such waveUke be- 
havior. If these solutions actually approach a final steady state, then the discus- 
sion in the previous paragraph applies to the actual solutions rather than their 
time averages since in this case the time average of the steady state is identical 
to the steady state. 

Section m. of this paper is devoted to the large time behavior of the first 
local approximation to the quasi-linear diabatic solutions. With some further 
assumptions on the properties of the magnetostatic turbulence which are made 
in this section and are further discussed in Appendix B, we are able to construct 
an H-theorem for the first local approximate solutions. From the H-theorem 
we find that these solutions must asymptote, with increasing time, to an iso- 
tropic steady state with a first derivative with respect to /x which is zero almost 
everywhere. The H-theorem does not require the first derivative in this final 
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steady state to be zero at = 0, and it actually allows the first derivative to be 
non-zero at other undetermined ^t-values as well. We continue by limiting our 
consideration to those solutions which conserve probability and which remain 
integrable in ix for all r. We also assume that the magnetostatic turbulence is 
such that the zero’th and first moments of the Laplace transformed kernel in 
the quasi- linear diabatic equation are continuous functions of jx with no zeroes 
in the first moment as a function of Then, we are able to show that the first 
derivative with respect to jx of the first local approximate solutions must be a 
bounded function of r pointwlse in /a. Thus, these solutions must be continuous 
in the closed /A-interval and bounded functions of p. pointwise in r. This infor- 
mation allows us to prove that the first derivative with respect to (x of the first 
local approximate solutions must be a continuous function of /a for any finite 
time, and when the time goes to infinity this first derivative must remain con- 
tinuous in jx except at ja = 0 where It is discontinuous and talces on the large 
value which we found the time average of the first derivative of the quasi-linear 
diabatic solutions must approach. Because we found from the H-theorem that 
the first derivative must asymptote, with Increasing time, to zero almost every- 
where in /A, and since here we find that for fx 0 the first derivative must re- 
main continuous in (x even as the time goes to infinity, we can finally conclude 
that the first derivative must actually asymptote to zero everj^here in ix except 
at ^ = 0. Thus, we find that the first local approximation to the quasi-linear 
diabatic solutions must asymptote, with increasing time, to exactly the same 
state which the time average of the quasi-linear diabatic solutions must also 
approach. 
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Given the very good agreement between the quasi- linear diabatic solutions 
and their first local approximation which is evident from our numerical results 
for short times, and given the mannei of agreement discussed in the preceding 
paragraph as the time asymptotes to infinity, we feel that it is reasonable to' 
assume that the first local approximation remains an excellent approximation 
for all time greater than zero. 

In section IV. of this paper we discuss a general consequence of the local 
approximation method which has broad implications to statistical mechanical 
studies of the approach to equilibrium. We restrict our attention to those solu- 
tions of the general integro-differential equation given in paper II. by equa- 
tion n.l and to those local approximations which remain bounded functions of 
time. We generalize our conception of a local approximation in this section to 
include any of the approximations which result when the Laplace transform of 
the kernel in equation II. 1 is replaced by any one of its Fade' approximants.'’ 
Then we show that the exact solution of equation n.l and all of the local approxi- 
mations which are generated by Fade' approximants that do not have a reduced 
range compared to the exact kernel transform have time averages {as defined 
through equation III.l) which must asymptote with increasing time to a single 
unique final state . As a corollary to this result we show that the failure of the 
adiabatic approximation can be understood on the grounds that it is in fact gen- 
erated by a Fade’ approximant which does have a reduced range. 

II. THE QUASI- LINEAR DIABATIC SOLUTIONS 

In this section we consider the final behavior of the quasi-linear diabatic 
solutions of equation 1.20 (paper I.). We introduce the following type of time 
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averaged probability distribution function, 


f>T • (:;r4^) I d.f(M,s) 


III.l 


and show that <f must asymptote to a final isotropic state (i.e., independent 
of /x) as T In this expression T is any finite time which, however, can be 

made as large as necessary so long as t is larger. Thus, equation in.l gives 
a time average of the final evolution of the quasi-linear diabatic solutions. We 
also consider the first derivative of r) with respect to /i . We introduce 
another time average. 


< 


f'>n 



df(^, s) 


III.2 


and show that at = 0, < f does not asymptote to zero as might be expected 
in view of the approach to isotropy of < f >j, but instead asymptotes to a large 
but finite value given by. 


< f '> 

* n 


1/2 - N(0, 0) 


tT 


III.3 


We also show that when /i. / 0, <f'>^~ Oasr - «. 

These results hold for "normalizable" quasi-linear diabatic solutions in magne- 
tostatic turbulence which has a "simple" power spectrum associated with it. 

Both of these terms are defined immediately below. 

In our numerical evaluations of the solutions of equation 1.20 from paper I. 
(the quasi-linear diabatic equation), we have never found evidence for negative 
values of f (/z, r) for any /j, and t. Nevertheless, equation 1.20 may contain solu- 
tions which do become negative. We limit our discussion here to those solutions 
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which are bounded from below so that a sufficiently large time independent 
isotropic component of the probability distribution function can be added on to 
the quasi-linear diabatic solution to make the sum non-negative, and thereby 
qualify it as a probability distribution function. Thus, having made the restric- 
tion to quasi-linear diabatic solutions which are bounded from below, it is no 
further restriction to assume that these solutions are non-negative. We further 
restrict our consideration to solutions which conserve probability. From equa- 
tion 1.24 of paper I., we see that we expect N (± 1, r) = constant so long as the 
integral on the right side of that equation remains bounded. Due to the finite 
range of the correlation function at /z = ±1 (see equation 1.21), we expect the 
integral to remain bounded unless 

t) _ B f (/z, r) 

B ^ /z 

becomes unbounded. We assume a bounded derivative at /z = ±1. Then N(±l, r) 
is a ccistant in time. Since equation 1.24 is homogeneous we further set N (1, r) 
= 1 with no loss of generality. We call solutions which satisfy the above re- 
quirements, "normalizable”. Notice that for normalizable solutions, 

0 < N(;CZ, T) < 1 III.4 

because f (mi > 0. In the last part of this section we will further limit our 
discussion to those solutions which have a bounded first derivative with respect 
to /z as a function of r in order to deduce the behavior of this derivative when 
Iji / 0 and T - 80 . 
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From equation 1.38 of paper I. we find, 


K(^, 0) = ^o(^) 



m.5 


in which P(oj) is the power spectral density, or more simply, the power spec- 
trum. As in paper I., we will assume a "simple" power spectrum which is non- 
negative and which does not contain any zeroes in the frequency range, 1 £ a. < «. 
In addition we will assume that 



(o) -• 3o) 


III.6 


so that the total power in the random field is finite. We note here that in our 
dimensionless notation, co = 1 corresponds to that wavelength in the magneto- 
static turbulence which is equal to the correlation length. In the next section of 
this paper, where the first local approximation to the quasi- linear diabatic solu- 
tions is discussed, another special requirement will be imposed on the turbu- 
lence; the simple power spectrum is sufficient here. 

The remainder of this section will be subdivided into four parts. In the 
first, we will develop certain properties of the Laplace modes and their asso- 
ciated eigenvalues which are necessary for the remainder of this section. In 
the second part we will prove that f asymptotes to isotropy, and in the third 
we will prove equation III. 3 nevertheless. In the fourth part of this section we 
will prove that <f '>^ ~ 0 as r - 0 when ^ / 0. 
a). The Laplace Modes 

We reintroduce the Laplace mode expansion of the probability distribution 
function and prove several useful properties of the modes and their associated 
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eigenvalues. On integrating equation 1.64 of paper I. in (m from -1 to ju, and 
then rearranging the result slightly, we find, 




^ p) 




0m (A^'- P) 


III. 7 


in which 


$ (At. p) = (1 - p) 


III.8 


where K(yu., p) is defined through, 

K(/i, p) = r dr R(/i r) cost 

Jo 


III.9 


(To obtain equation Ilf.? we have divided by (p). Since (p) =0, we exclude 
m = 0 in the following. We also point out that the divergence of the right side 
of equation III.7 at yu. = ±1 i:: only apparent; the integral also goes to zero at 
these points so that the ratio remains finite.) The quantity, R(0> is the two- 
point correlation function associated with the magnetostatic turbulence (see 
equation 1.19) and is a function of the spatial separation, C- 

Having adopted the simple power spectrum, we find that when /j, 0, 


p) Pi 


(/^ ^ 0 ) 


III. 10 


and therefore that, 


^n.(P) 


^0m(M. P) 




Pi 


(M 0) 


III.ll 


On the other hand, when ya = 0 we have. 
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K(0, p) 


p 


m.i2 


1 + p2 

(we have normalized the correlation function so that R (0) = 1) and, in contrast 
to equation ni.lO 


Thus, 


P 

®(0. p) 


Pi 


(M = 0) 


III. 13 



3 0,„<O>P) 

. — 

Pi 



UI.14 


In paper I. (see discussion surrounding equation 1,69) we have shown that the 
{j. 1 , 0) can be constructed through even and odd combinations of the eigen- 
fimctions which form a complete orthogonal set on the half- domain , - 1 1 Ai- 5 O' 
In this half-domain, and starting at m = 0, each pair of (ju., 0) (e.g., m = 0, 1; 
m = 2, 3; etc.) are identical to each other. Thus, from the orthogonality of 
these eigenfiuictions on the half-domain we can conclude that the integral on the 
right side of equation 111,14 is non-zero for m = 0 and 1 only. Therefore, 


P 

K(p') 


P) 


-1 


PI 


0 


(m = 1) 

(m > 1) 


We have further seen in prper I. that, 


P) 

■' 'V 

Pi 


cc 


III. 15 


HI. 16 


Thus, even though (p) 0 (p - 0), we still have. 
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p 


0 


m.i 7 


^i(p) Pi 

For m > 1 we found in paper I. that ^j(p) ~ (p - 0) with all > 0. 
Therefore, we can make the general statement, 


(P) Pi 


(m > 0) 


III. 18 


Finally, from equation III.7 we find. 



(ax, p) ^ 

V ^m(P) ) 



’/'h.Cm'.p) 


IU.19 


and therefore. 


P \ P) / $(m,P) \ ^ 

\(P)/ \'^„(P)/Pl 


m > 0 


ni.20 


With this information, we can now proceed to a discussion of the final states of 
the time averaged quantities introduced phove. We note in passing, that we have 
only assumed the simple power spectrum to obtain the results presented so far 
in this section, 
b). Isotropy 

From equations 1.65 through 1.67 of paper I., we find that the Laplace mode 
expansion of pf(/Lt, p), where f(/i, p) is the Laplace transform of the probability 
distribution function, is given by. 


pf (m. p) 




f P) f (^'>0) 

E.^(P) 

* t.(pj] 



III. 21 
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where. 


fo(0) ■ \ fW'f(M'.O) = j 


UI.22 


the last equality holding for normalizable solutions. 

By construction, and through the properties of the Laplace modes which are 
solutions of the Sturm- Liouville problem (equation 1.64, paper I.), we are guar- 
anteed that this expansion converges absolutely and uniformly in /j. for any p >0. 
Furthermore, by our construction of this same expansion for f (^t, p) at p = 0 in 
terms of the adiabatic Laplace modes we are also assured that this expansion 
converges absolutely and vmiformly in jx for p = 0 as well. From our discus- 
sion in section a), concerning the properties of the eigenvalues as p -♦ 0 we 
conclude, 


pf(/i, p) ~ fo(0) 

Pi 


m.23 


Because the normalizable solutions are non-negative we can apply the Tauberian 
theorem, theorem 4.3, Chapter V., of Widder, which is a special case of a 
theorem provided by Karamata, as well as a special case of Wiener's general 
Tauberian theorem, and which allows us to conclude from equation in. 23 that. 


u: 


ds f (/X, s) fo(0) 


m.24 


Tt 


From equation in.l we can obtain, 
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<f>. 


Jo 


ds [ (yu, S) - 


^ I ds f (/X, s) 
Jo 


m.25 




f 

J n 


Therefore, 


<f>T ~ fo(0) ni.26 

tt 

for any T less than infinity. Thus, we conclude that <f must asymptote to an 
isotropic state, with increasing r, which is actually the isotropic component of 
the initial probability distribution fimction at t = 0. 
c). The First Derivative Qx = 0) 

We introduce the notation, 

Q(r) - f ds Q(s) m.27 


for any quantity, Q, and note that if the integral on the right side of this expres- 
sion is bounded for any finite r, then Q(t) ~ <Q^t» “* “)• 

We also reintroduce the functions g(/x, t) and h(/x, r) which were originally 
defined by equations 1.47 and 1.48 in paper I. We adopt these functions here with 
the exponential correlation function replaced by an arbitrary functional depen- 
dence on spatial separation. With g(yLt,r), the quasi-linear diabatic equation for 
N(/i,r), equation 1.24 of paper I., can be written. 


3 N 
3 T 




III.27 


On differentiating t) and h(;u, t) with respect to time we find. 


14 



- 

d T 
and, 


(1 -;x2) 


B 


f 

^'0 


, . 2 x r »R(a‘(’- - M) , ,, 3^N(\) 

h + (1 - /j. ) I d\ cos (t - \) 

o T 3 ^2 

m .28 


Bh 

B T 


B R(/x(t - \)) 


B2N(\) 


^ = g + (1 - M^) d\ ^ ^ sin (T - \) ' ' ni.29 

Jo 


B T 


B 


At ^ = 0, these equations reduce to 


BN(0, r) _ 2 


B T 


= 7? g(0. '^) 


ni.30 


Bg(0. T) _ B^NCO.r) 


B T 


Bu2 


- h(0, T) 


in.3i 


and 


Bh(O.T) 

= g(0,r) 


m.32 


because, 


3r(;u(t - X)) 


= mR'(0 


> 1-0 

£ =/ l ( T - A ) 


III.33 


for bounded R' (0) which we will assume. 

On applying the integral operation given by equation m.27 to equations UI.30 
through in.32 we obtain, 


N(0, r) - N(0, 0) 

^ -—Am 


7? g(0. 7-) 


III.34 


g(0, r) _ B2N(0 ,t) _ 

— - - h(0,r) 

^ Bu2 


111.35 
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and, 


h(0, T) 

T 


i(o.T) 


III.36 


Through the substitution of equation in.36 into equation III.34 we find, 


h(O.T) 


N(0,r) - N(0, 0) 


m.37 


from which we conclude that h (0, r) is a bounded function of time in view of 
equation in. 4. From equation ni.37 we determine h(0, t) and substitute the 
result into equation ni.35 to obtain, 


g(O.-r) 

T 


N(0, T) / N(0, t) - N(0, 0) ' 
3^2 ^2 


IIL38 


Our evaluation of the asymptotic behavior of this expression depends on the be- 
havior of g(0, r) as T no. 

From equation in.32 we see that, 

h(0, r) = r ds g(0, s) m.39 

Jo 

Since h(0, r) is a bounded function of time we conclude that g(0, r) must also be 
a bounded fimction of time except that it could contain integrable singularities. 
However, it is possible to show that if we assume the existence of an integrable 
singularity in the time dependence of g(0,r) we obtain a contradiction with the 
boundedness of h(0, r) through its definition which is. 


h(0,T) 




dX. sin (t - K) 


B^ N(0, \) 
B 


m.40 
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If g(0, r) contains an integrable singularity, then from equation III.39, 
h(0,T) remains bounded at the position of that singularity. However, if g(0,'^) 
contains an integrable singularity then g(0, t) contains a non- integrable singu- 
larity at the same position in r. Since h(0, r) is a bounded function of time we 
see from equation 111.31 that in this situation N( 0 ,r)/ 3^2 m^gt also contain 
a non- integrable singularity in time at the position of the singularity in g(0, r). 
But, from equation III. 40 we see that then h (0, r) must become unbounded at the 
position of the singularity leading to a contradiction with the assumed integra- 
bility of the singularity in equation III.39. Thus, we conclude that g(0, r) must in 
fact be a bounded function in time. Therefore, 


TT 


1/2 - N(0, 0) 


in.4i 


follows from equation HI. 38 as t - 
d). The First Derivative (/z 0) 

The Laplace transform of equation 1.24 (paper I.) is. 


pN(/x, p) - N(/x, 0) = 


3 f (m. p) 


Therefore, 


3f(Ai, p) p[pN(/i, p) - N(/i, 0)] 


3/x 


J (/X, p) 


m.42 


m.43 


Division by ID (/x, p) is allowable since at = ± 1, p N(/x, p) - N(yx, 0) = 0, and 
the right side of this equation remains bounded. In addition we exclude /x = 0 
from consideration where f (0, p) has a zero when p = 0. From equation III. 10 
we see that we expect the right side of this equation to ~0 as p "* 0 unless the 
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quantity on the right side in the square bracket becomes unbounded as p - 0. 

In the following we first show that this quantity in fact remains bounded as p 0. 
We then use this result to deduce the behavior of < f as r - « when ^ 0. 

Since the infinite expansion given by equation III.21 is known to converge 
uniformly in yu, we can integrate that expansion term by term to find pN(/z, p). 
We find, 


pN(/i, p) = (/i + 1) f(,(0) + 


j P) f(M'. 0) 


/ 


E„“(P) 


m- 1 


^m(P) 


where, 


P) * ) clAi'0„ (m', P) 


m.44 


ni.45 


But, from equation III.7, 


Therefore, 


^m(^'P) - 


/ £ (M. P) \ 
\ ^m(p) / 


pN(/i, p) - (^ + 1) 1^(0) 


111.46 


oc 



m- 1 


<\(p)/ xKip)/ J-i ' 


EJ'(P) 


+ 


K(p)\ 


m.47 
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and, from equation III.20, we see that, 


pN(M. P) - (Ai+l)fo(0) III.48 

Pi 

which is certainly bounded. Since N(fx, 0) is also bounded we conclude that 


a f (m. p) 


~ 0 
Pi 


III. 19 


We now assume that B f (^, r) / B;x is a bounded function of time. Therefore, we 
can apply Wiener's generalized Tauberian theorem to Laplace transforms, lo 
obtain (see Appendix A), 


i r ds 

*'o 


3 f (m. s) 

OfM 


m.5o 


tt 


and then, we can conclude. 


<f'>T ~ 0 (fjL ^ 0) ni.51 

TT 

Before proceeding, we note that as a corollary of the above results we can con- 
clude that if the quasi-linear diabatic probability distribution function actually 
approaches a final steady state given by. 


and 


then, 


f(M, '^) ~ S(/i) 

tt 


B f (ll, t) 

,, I 

B/i 




S(/x) = fo(0) 


III.52 


III.53 


III.54 
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and, 


l/2-N(0,0) ^ __ 

- 0 ) 

V 

S'(fji) = m.55 

0 (At 0) 

otherwise, as we showed above, this final steady state is the time average of the 
final evolution of the quasi-linear diatatic solutions. One of the main conclu- 
sions of the next section will be a proof that the first local approximation to the 
quasi-linear diabatic solutions must approach this final steady state as t - 
(no time averaging is necessary). 

m. THE FIRST LOCAL APPROXIMATION TO THE 

QUASI-LINEAR DIABATIC SOLUTIONS 

In this section we consider the first local approximation to the quasi-linear 
diabatic solutions whose large time properties were discussed in the preceding 
section. The central issue here is an H-theorem which we will construct for 
the first local approximation and which we will use to show that in very general 
circumstances the approximate solutions must asymptote to the steady state 
characterized by S(/i) and S' (jj.) above. Thus, in this section w'e show that the 
approximate and exact solutions must asymptote to exactly the same final steady 
state if the exact solutions have a final steady state. If the exact solutions do not 
have a final steady state, then what we show in this section is that the first local 
approximations and the time average of the exact solutions must still asymptote 
to that final state given by S(u.) and S' (m). 

We divide this section into four parts also. In the first part we construct 
the H-theorem and from it deduce that f ^(/x, t) isqtropy, (t - «>) and f j' (/x,r) 

0 (r - oo) almost everywhere in ^ (I I (/x, r) = 3fj (/x, r)/ 3^). The amplitude 
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of the isotropic final state is not given by the H-theorem. However, we restrict 
our consideration to solutions for which f^' (± 1, t) and B fj' (± 1, t)/Bt are 
bounded, so that total probability is conserved, and then the amplitude is given 
by S{fj.) = fg (0). The H-theorem does not determine the asymptotic behavior of 
fj' (/i, r) at those points where f^' 0 (t - »), it does not tell where those 

points are, and it also tells nothing of the asymptotic behavior of f j' (0, t). The 
remaining three parts of this section are directed to a fuller understanding of 
these issues. 

In the second part of this section we restrict our attention to those solutions 
for which f ^(/a, t) remains integrable in so that N j r) exists and is finite 
for all (/i, r). Since this property is an essential ingredient in the definition of 
a probability distribution function we feel that this restriction is reasonable. 

For these solutions we show that f (;x, r) is a boimded function of time point- 
wise in y.. From this result it follows that both f^ (yu, t) and Nj (/x, r) are con- 
tinuous functions in the closed interval - 1 < M < 1, and therefore that f^ {y., t) 
is a bounded function of y. for any t. Thus, these solutions are normalizable in 
the sense of section II. In the third part of this section we show that because 
these solutions are normalizable, 


fi'(0,T) 


1/2 - N(0, 0) 


TT 


HL56 


and in the fourth part we show that fj' (y, r) is continuous in y. for any finite r, 
is continuous in y at t = « and /x ^ 0, and is discontinuous in /x at r = « and 
y = 0. To prove this result we must make the further assumption that both 
iTg (;Li) and (y) are continuous fimctions of y. The significance of this result 
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») for all ^ 0. Therefore, 


is that we can then conclude that f^' (a^, t) ~ 0 (r - 
we conclude that (^t, t) ~ S' (/x) (t - «). 
a). The H-theorem 

In order to construct the H-theorem we require, 

iHli(At) = r dr R(/zt) t cos r < 0 III.57 

with lUj (/u = ± 1) bounded. This requirement is met by the exponential correla- 
tion function used in paper I. and n. to obtain the numerical results given there. 
However, we have not been able to devise any simple criterion for the correla- 
tion function or its associated power spectrum in order that equation III.57 be 
satisfied. For the remainder of this section we do assume that this condition on 
the magnetostatic turbulence is met, and in Appendix B we list a few alternative 
expressions for equation ni.57 none of which, unfortunately, are particularly 
transparent. We also continue to assume a simple power spectrum as defined 
in the previous section. 

Equation 11.62 (paper U.) generates the first local approximation to the 
quasi-linear diabatic probability distribution function. We reproduce this equa- 
tion here for f^ (ai, r) instead of Nj (al, t): 


B r 



(m) 


111 

BAi 



IIL58 


subject to the reinitialized condition which is given in terms of the initial condi- 
tion by. 


, B ^ Bfj(Ai, 0^) 

f,(M, 0") + ^ = 0) 

d At d At 


m.59 
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We further introduce the quantity, H(r), through, 


H(r) 


1 

2 









III.60 


and note that through its definition and the assumed properties of H(r) >0. 

Through integration by parts we find that the time derivative of H (t) is given by. 



- 3 

- 3 - ^ fi 

H(r) = 

• 

Tr 

.1 

a fj. ^ a fj. 


III.61 


I 


= - -r.2 


dM ^0 (m) 



J-l 

and thus, H(t) < 0. It is impossible for H(r) to remain greater than zero 
unless H(t) ~ 0(t - <»). Thus, we conclude that fj' (/j., r) ~ 0(r - *>) almost 
everywhere in /j.. Since IqIO) = 0, it is clear that f^' (0, t) does not have to 
asymptote to zero. It is also possible that f^' (^, r) does not asymptote to zero 
at other points in ^j. so long as these points do not contribute to the integral in 
equation 111,61. Nevertheless, we must have, 

= f,(Mo.r)+ Td/x' - fi^.r) m.62 

Ju 

where is any point in /z at which f j' r) does asymptote to zero. Thus, we 

conclude that f ^ {ji, t) must asymptote to an isotropic function of t. But, it is 

* 

also easy to see that this function must be constant in r . Because H ( r) ~ 0 (r - *), 
we must have H (t) ~ (r - «) where H ^ is indeed a constant number which is 
independent of t. But, from equation UI.60, we see that H(7 ) ~ f ^ r) (t - «). 
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Therefore, “* “)• conclude that f ^ (/x, t) must asymptote 

to an isotropic steady state as r -* e». We further restrict our attention to those 
solutions for which f^' (± 1, t) and B f (± 1, t)/Bt are bounded functions of r. 
Then we see from equation n.62 (paper II.) that (± 1, r) is constant in time; 
i.e., these solutions conserve total probability. In this case it is clear that 
fi(M. t) ~ fo (0) (r ^ »). 

hi the remainder of this section we determine the sufficient conditions to 
ensure that f^' {jx, r) ~ 0 (r oo) everywhere in /x except at /x = 0. 
b). fj’ {jj., t) Is Boimded Pointwise in fj. 

We have already assumed that i ^ (± 1, r) is bounded in r . Here, we show 
that f j (jjL, r) is bounded pointwise in ^ in the open interval, - 1 < yu < 1. To 
proceed, we define, 

q(yLX,T) S _ <^2 jjj gg 


and, 


A(yu) - - 


®1 


> 0 


(-1 < /X < 1) 


m.64 


Here, we are assuming that JHj(m) has no zeroes in the open interval. Then, 
equation 11.62 (paper II.) can be rewritten for positive r as. 


^ q(M. t ) 
?r 


+ q(M. '^) - 


B Ni(/x, T) 
B T 


m.65 


and integrated. 
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q(M. 


m.66 


- e 


- Ay-" 


q(^, 0+) + f ( 

Jo 


dK e‘‘^ M A 




With the use of integration by parts and reiniiialization we find, 

q(/x, r) = N, (/X, t) - e" N(/u., 0) - A(/i) e~ f d\ ^ (/u, \) 

Jo 


m.67 


Thus, 


qC/i. '^)1 1 iNjCAi, ' t)| + le" N((U., 0)1 


III.68 


A(/i) r d\ eA(M)A 

^0 


We have already assumed that the initial condition is normalizable and that 
Nj (/X, r) is imiformly bounded in ^ and r ; i.e., | (/u> t)| < where is 

finite. Therefore, 


lq(A^. ^)l 1 Ml + 1 + Mj 1 1 - 

< 2Mj + 1 


III. 69 


From equation III.63 we see then that f j' (/x, r) is bounded on the open interval 
in /j.. We can further conclude that both f^ (/j., r) and N j(/x, r) are continuous on 
the closed interval, - 1 £ /x < 1. Therefore, both of these quantities are uni- 
formly boimded in /x. We are thus able to assume that these solutions are nor- 
malizable without loss of generality, and that. 
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0 < T) < 1 (-1 < 1 1) 

for any r. 
c). fi» (0 ,t) 

Since IDq (0) = 0, we see from equation 11.62 (paper II.) that, 

— Ni(0, T) + 772 !D,( 0 ) — =0 (r > 0) III.71 


Therefore, the quantity in brackets in this equation is a constant in time which, 
by integrating equation in. 59 in /x, we can find is N (0, 0). Thus, 

B?i(0, r) 

N,(0, r )-772 = N(0, 0) m.72 

dfi 

since IDj (0) = - 1. Given that these solutions are normalizable, we have 
Ni(0, T) ~ i/2(- - Therefore, 



whi'^^' is exactly the same value that this derivative must talce on in a time 
averaged sense for the exact quasi- linear diabatic solutions, 
d). / 0 

It is relatively easy to see from equation III.67 that q^u., t), and therefore 
f j {fjL, t), is a continuous fimction of jj. for any finite r, and when t = »> that 
q{(i,r) remains continuous in /x for jj. ^ 0. We examine | q(/Xp r) - q(/XQ, r) | 
where and llq are two different values of fi, and show that this difference can 
be made as small as we like if yxj - fj.^ is made small enough e.\cept when 
yxg = 0 and T = oo (we hold /Xq fixed and let - /Xq ). 
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We define, 


e = Ni (/^i, t) - t) 


III.74 


and, 


S - A(^^)-A(/Xo) 

We have already shown that (/x, t) is a continuous function of fi, and we have 
assumed that A (/x) is also continuous. Therefore, both e and S can be made as 
small as we like for small enough In general e should be considered 

a function of time, however, in light of equation IIL70 a bound on this difference 
can be fotmd which is uniform in t. Thus, 

lq(/xj, r) - q(/XQ, r)l < | (^ii, r) - (/x^, r)l 

+ (e-S^ - 1) N(/Xo, 0) + e 

III.76 

+ A(mo) (e-®^ - 1) r - \) 

•>0 

+ |e| A(;xo) +|S| f d^ e“ ^ N (/x^, r - \) 

Jq -’o 


and then, with the use of equation III.70, 
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iC^i, T) - q(/Xjj, r)j < |g| + (1 - g-Sr^ 


+ Id 


^(Mo) 


!l - g“ (‘^^(Mo) + S)r 

V" ^(Mo) + S 


j - (l - e-'^<"»>") 


+ Ul (1 - + |s| \ 

\ / 


III. 77 


It IS clear that for finite r , this expression can be made as small as is neces- 
sary by making e and S small enough When r - « we consider the two cases, 
Mo 0 and ^0 = 0, separately. If 0, then (we assume ^ 0 also), 


lq(Mi,«) - q(^o,«o)| < 2 |e| + 2 


^(Mt) 


m.78 


which can be made as small as we like. However, if =0. then, 


lq(Mi- *) - q(0, «)1 < 3 + I e| 


III. 79 


In this case we cannot make the difference as small as we like by reducing j 
and 5 . Thus, we conclude that q , t) is continuous in for any finite t , and 
that q(^, r) is continuous in for r = „ when m 0; q<,,, r) is not continuous 

at M 0 and r - Since f,, (^) has been assumed continuous in a , these same 
conclusions hold for f/ (^, t). 

Since in a), we have learned from the H-theorem that f/ (^ , r) is zero 
almost every^vhere at r = », we can now conclude from the continuity of fj' (^, r) 
that in fact fj' (/x, r) ~ 0 (r - «) for all ^ except yx = 0. We have seen in c). that 
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f i’ (Ml '^) does not asymptote to zero at /x =0; this is consistent with equa- 
tion m.79. 

We conclude, finally, that under the conditions outlined in this section, the 
first local approximation to the quasi-linear diabatic solutions must asymptote 
to the final steady state characterized by S(^) and S' (;i) in the previous section. 
As peculiar as this final state is, we see that the first local approximation yields 
exactly the time averaged behavior of the quasi-linear diabatic solutions as 

T -* oo , 

In Figure 1, we have plotted f i' (m, '^) at various r for the broad parallel 
beam injection considered earlier in papers I. and II. in order to give some in- 
dication of how' the final state, S(u) and S' (/x), develops; figure 20 of paper n. is 
also useful for this purpose. As r increases, we can see from figure 20 of 
paper II. that f^ (/x, r) tends to the isotropic state, S{^) = f^ (0), but with a nar- 
row region near = 0 in which a steep gradient persists. From Figure 1 of 
this paper, we see that the width of this region decreases with increasing r 
while the amplitude of the gradient at = 0 increases. We suppose that as 
T - the amplitude of the gradient saturates at the value given by equation III.73, 
the width of the peak in the gradient shown in Figure 1 goes to zero, and f j (^x, r) 
reaches isotropy. The amplitude of the peak in the gradient is plotted in Figure 2 
as a fxmction of From this figure, the approach of this peak value to the 
maximum S' (0) seems reasonable. We also point out that the approach to S' (/x) 

= 0 (/XX ^ 0) is consistent with Figure 1 for jx ^ Q. 
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Figure 1. The first derivative with respect to ^ of the first local ap- 
proximation to the broad parallel beam injection solution computed 
earlier in papers I. and II. The derivative is plotted as a fimction oi/j. 
for several values of time in units of gyro-periods. The irregularity 
of these curves is due to the presence of numerical noise which has 
been amplified by the differentiation process. 
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f'(O.oo) 


BROAD, PARALLEL 
BEAM INJECTION 


Figure 2. The first derivative with respect to /x, evaluated at m = 0, of the 
first local approximation to the broad parallel beam injection solution computed 
earlier in papers I. and II. 




IV. THE GENERAL N'TH LOCAL APPROXIMATION 


In paper II we constructed the local approximation method as it applies to 
the integrodifferential equation defined by equation II. 1. The Laplace transform 
of that equation is, 

pf(p)-f(0) = eK(p)f(p) m.80 

and, the Ljqjlace transform of the equation which generates the n'th local ap- 
proximation to its solution is, 

Pf„(p)-f(0) = ^K„(p)f„(p) m.81 

where K„ (p) is the n'th convergent (see paper II, section II) to K(p). A formal 
solution of either of these equations can be written as, 

f„(p) = [p - e K„(p)]"‘ f(0) ni.82 

The exact solution of equation HI. 80 is obtained in some exceptional cases for 
all n > N (N finite), but in general can be obtained by letting n - »=, Thus, 
equation III. 82 contains all of the local approximations, as well as the exact so- 
lution, within the entire set of n-values. 

We have retained the notion, and associated notation, here that the local 
approximations are generated by those entries in the Pade' table (see paper II.) 
which are near the diagonal and are the convergents to K (p). However, the fol- 
lowing discussion holds true for those approximations to the solution of equa- 
tion II. 1 which are generated by a^ of the Pade' approximants. Thus, for the 
discussion in this section we can generalize the definition of the local approxi- 
mations to include all approximations which are generated by any Pade' table entry. 
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From equation III.82 we see that generally, 




-1 

“ / \ * 





pf„(p) = 

)\ 

f(0) ; 
Pi 



f(0) 


ni.83 


for any n. Thus, through the application of the Tauberian theorems which we 
discussed in the previous two sections we can come to the very important con- 
clusion that the time averages (as defined through equation ni.l) of all of the 
local approximations, including the exact solution, must approach a single final 
state given by. 



m.84 


There are exceptions to this very general statement; the failure of the adiabatic 
approximation is due to one of them. In the following we will discuss these ex- 
ceptions briefly. 

a). The Exactly Soluble Model Equation of Paper II 

The exactly soluble model equation has a C-number kernel whose moments 
are given by equation 11.48. In this case, from equation III. 84, we expect that 
the time average of all of the local approximations, including the exact solution, 
to asymptote to zero as t -* ». From figures 2 through 5 of paper II we see 
that this expected behavior is obtained except in some special situations. First, 
from figure 4 we see that the third local approximation does not fit in to the 
general pattern predicted above. Indeed, as was discussed in p^er n, 13 (t) 

~ - *1 (r - «) for 0 < M ~ 0.6. Since this local approximation is not bounded 
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from below, we simply cannot apply the Tauberian theorems which lead to equa- 
tion ni.84. Thus, there is no contradiction here; we cannot apply equation III.84 
to this case in the first place. 

The remaining exceptions to the general statement made in the previous 
paragraph occur in figure 5 at yit = 0 in the exactly soluble model equation. In 
this case we see that the time averages of the first, second, and third local 
approximations agree exactly with the time average of the exact solution, but, 
none of these time averages are zero. In addition, the time average of the 
zero’th local approximation (the adiabatic approximation) obviously does not 
agree with any of the others. This peculiar behavior occurs because Mq (0) =0, 
and Mj (0) = - 1 ^ 0. The time average given by equation III. 84 must be re- 
placed by, 

^ f(0) m.85 

in this situation. Since M j (0) = - 1 and 6 = 1 in the example being discussed 
here, we see that application of equation III.85 does lead to the prediction of the 
correct time averages (= 0.5) for the first through fourth local approximations. 
On the other hand, we see that the zero'th local approximation must fail to pro- 
duce the correct time average at /Lt = 0 because the dominant term in the con- 
tinued fraction expansion of K(p), namely Mi (0), is apriori neglected in this case. 
Thus, we see that equation III.84 does yield the asymptotic state of the time 
averages of all the local approximations so long as 0. If = 0, then 

other considerations must be made which lead to replacements for equation IU.84 
. . . for example, when Mj / 0 then equation III.85 follows . . . and the adiabatic 
approximation must be abandoned. 
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b). The Quasi- Linear Diabatic Equation 

We have seen in paper I (equations 1.65 through 1.67) that the Laplace trans- 
form of the exact solution of the quasi-linear diabatic equation can be written, 


p) - X! (p) '/'m P) 

m=0 

where the >/'„(At, p) are the Laplace modes and where, 

(/X, p) 7(/i, 0) 


III.86 


Pfm(p) = 




m.87 


1+7)2 


\ (p)'^ 


In addition, we found there that \q(p) = 0, and that \j(0) = 0, but, from part a), 
of section II in this paper we foimd that p (p) ~ 0 (p - 0) for all m > 0 in- 
cluding m = 1. Therefore, pf„,(p) ~ 0 (p - 0) for all m > 0 and 


1 - 

pfo(P) ~ 0 dfz f (fj., 0) 

Pi ^ 


III.88 


Thus, through the use of the Tauberian theorems we conclude (as we have 
earlier in this paper) that the time average of the quasi-linear diabatic solution 
must asymptote to the final isotropic state, 


- [' 
2 J.i 


djj. f (fi, 0) 


m.89 


Equation III.89 gives the projection of the initial state onto that Laplace 
mode which is the only one outside the range of K(p). The zero'th Laplace mode 
is given by » P) ~ 1* the projection of the initial state onto » P) 
given by equation in.89. We say that p) is outside the range of K(p) because. 
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K(p) P) = '^o(m. P) = 0 IU.90 

It is that projection of the initial state which is outside the range of K (p) which 
survives as t -« «> and becomes the final state which the time average of the 
solution must approach. The first Laplace mode is not outside of the range of 
k(p) since, 

K(P) v^iCm-P) = -^i(p) '/'i(M.p) ^0 (p > 0) m.9i 

But, in the zero'th local or adiabatic approximation we replace this relationship 
with, 

Kq(p) 0) = -\i(0) 0) = 0 m.92 

with the hope that 0i(m. 0) and (0) are reasonable approximations to '/'liM. p) 
and \j (p) when p is small. Unfortunately, in this attempt we have introduced a 
convergent to K(p) which, crudely speaking, "contains a zero" which the original 
K (p) did not. (This feature is analogous to that of the simple model equation 
discussed above in which (0) = 0.) Speaking more precisely, we have intro- 
duced a convergent to K (p) which has a reduced rai^e compared to K (p). This 
reduction in range then allows additional projections of the initial state outside 
the range of the kernel transform which survive as t - so and falsely modify 
the final state of the time averaged solution. More specifically, we have seen 
above that pfj (p) ~ 0 (p - 0) and therefore the projection of the initial state 
onto i/'i (^, 0) cannot contribute to the final state. However, if we replace \j(p) 
by its zero'th local approximation, (0) = 0, before taking the limit, p - 0, w'e 
find erroneously that. 
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pfl(p) 


III.93 


1 

0) f (//, 0) 

Pi 2 J_i 

and then the erroneous final state for tho time average of the solution contains a 
step in it at /X = 0 due to the presence of 0) (see figures 2 and 6 of paper I). 

At this point we are able to make the statement preceding equation III. 84 
more precise. We can say that the time averages of all those local approxima- 
tions, including the exact solution, which are bovmded functions of time must 
approach a single final state. This statement does not apply to those local 
approximations which are generated by a co nvergent, (p), which has a re- 
duced range compared to the range of the exact operator kernel transform, K(p). 
We feel that this general principle should be considered in future statistical 
mechanical studies of the approach to equilibrium. 

V. CONCLUSION 

In this series of three papers we have addressed two separate issues. In 
order to clarify these individual issues we present here a brief summary of our 
work. 

We begin our summary with the equation which Klimas and Sandri derived 
(see equation 1.8) from first principles, using the method of Kaufman, which 
governs the exact evolution of the ensemble averaged test particle probability 
distribution function in magnetostatic turbulence. This ''master equation" . . . 
an unfortunate name, since it has nothing to do with the well known master equa- 
tion formalism . . . follows directly from the Liouville equation through aver- 
aging over an ensemble of stochastic magnetic fields. In principle this master 
equation yields an exact solution; in practice, it must be simplified in order to 
extract any information from it at all. 
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Two steps must be taken in order to reduce the master equation to the 
Fokker- Planck equation. First, if either the magnetostatic turbulence is weak 
or if the energy of the test particle is high, a small parameter enters the master 
equation which allows its expansion. If a power series expansion of the master 
equation is constructed and then truncated beyond the leading non-trivial term, 
then the quasi-linear approximation results. The product of this truncation is 
an integro-differential equation in time which, in paper I., we named the quasi- 
linear diabatic equation (see equation 1.20 for a specific example in a simple 
model of the magnetostatic turbulence and also following gyro-phase averaging). 

In the second step, we apply the adiabatic approximation which reduces the 
quasi-linear diabatic equation to the quasi-linear adiabatic equation that is also 
known as the Fokker- Planck equation. In this step we have reduced the quasi- 
linear diabatic equation, which is integro-differential or non-local in time, to 
its Markovian approximation. The end result is a velocity space diffusion equa- 
tion which is first order partial differential in time and which governs the evo- 
lution of the probability distribution function in the spirit of the Markov chain; 
i.e., changes in the probability distribution function depend only on its instan- 
taneous state, and are independent of the manner in which this state was achieved. 

As we have discussed in paper I., solutions of the Fokker- Planck equation 
are so blatantly unreasonable that even in the absence of any exact solutions of 
the master equation we are able to conclude that something went wrong in our 
two step simplification which lead to the Fokker- Planck equation. The standard 
response to this failure has been various attempts, labeled non-linear, or orbit 
modification, or diffusing orbit, or partially averaged theories which are all 
attempts to renormalize the power series expansion of the master equation in 
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such a way that, on further applying the adiabatic approximation, the result is a 
velocity space diffusion equation much like the Fokker- Planck equation but with 
a diffusion coefficient which is considered more reasonable. These attempts 
are all based upon the apriori prejudice that a reasonable approximation to the 
unavailable solutions of the paster equation must be governed by a Markov chain 
of events as expressed through a velocity space diffusion equation similar to the 
Fokker- Planck equation. However, as we have just mentioned, exact solutions 
of the master equation are not available, and furthermore, to our knowledge, 
solutions of the quasi-linear diabatic equation, which follows after step one 
mentioned above, and before the imposition of the Markovian approximation, 
have also never been calculated prior to this work. Thus, the predictions of the 
non- Markovian theory have been imknown. In addition, as we showed in paper I., 
the imposition of step two, the adiabatic approximation which leads to the 
Markov chain, clearly leads to an incorrect description of the quasi-linear 
diabatic solutions. In fact, although it is only step two of the simplification 
procedure which leads to the Fokker- Planck equation which is clearly in error, 
it is only step one which has been modified in all previous attempts to correct 
the Fokker- Planck predictions. Apparently, the idea of a Markov chain is so 
imbedded that its abandonment has not been considered even in situations where 
it clearly introduces serious errors. 

Thus, we come to the two separate issues which we have considered in these 
papers. First, we ask, "Wdiat are the predictions of the non-Markovian quasi- 
linear diabatic equation, and do its solutions indicate that a renormalization of 
the power series expansion is indicated?" and second, regardless of the neces- 
sity or not of the renormalization, we ask, "Is there a modification, or 
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replacement, for the adiabatic approximation which then allows us to produce a 
successful approximation to the solutions of integro-differential equations like 
the quasi-linear diabatic equation?" After all, if modifications of the quasi- 
tinear diabatic solutions are necessary, then it makes sense to start with the 
quasi-linear diabatic solutions and see what direction these modifications must 
talce. In addition, since solutions of the quasi-linear diabatic equation are very 
difficult to obtain, even numer''‘'’’.ly, a simpler method for obtaining good approxi- 
mations to these solutiOiis cannot hurt in this endeavor. If this method happens 
to be powerful enough to be applicable more generally, then all the better. 

In paper n. of this series we introduced the new local approximation method. 
In that paper we showed that this approximation method can be applied to the 
fundamental equatii 's of statistical mechanics if they are cast in the integro- 
differential form of the quasi-linear diabatic equation (the master equation is 
included). In addition to -ihowing the wide applicability of this method we also 
demonstrated its application to a very simple exactly soluble problem and 
furthermore, to the quasi-linear diabatic equation which we constructed in 
paper I. In general we foimd this new approximation method to be remarkably 
successful. However, in the case of the quasi-linear diabatic equation, since 
both the exact and the approximate solutions could only be obtained through 
numerical integrations of the respective equations, the comparisons we have of 
these solutions extend only over a little more than two Larmor periods of 
elapsed time. Thus, in this paper we have presented an analytic study of the 
long time behavior of both the quasi-linear diabatic solutions and their local 
approximations. 
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We have found in this paper that the first local approximation to the quasi- 
linear diabatic solutions probably remains remarkably accurate for all time. 
More specifically, we have foimd that an H-theorem can be constructed for the 
first local approximation which, imder fairly general conditions, allows us to 
conclude that the first local approximate solutions must asymptote, with increas- 
ing time, to exactly that final state which the time average of the quasi-linear 
diabatic solutions must also approach. Given the remarkable agreement between 
these solutions during early times, and also this manner of agreement for very 
large times, we feel reasonably confident in the accuracy of the first local ap- 
proximation for all time. Thus, we feel that the second of the two questions 
posed above has been answered. The non-Markovian local approximation method 
does successfully replace the adiabatic approximation in applications to integro- 
differential equations like the quasi-linear diabatic equation. 

From the results of this paper, the first question posed above can now be 
answered as well. From paper I. (see discussion leading to equation 1 . 10 ) we 
can see that the power series expansion of the master equation which we dis- 
cussed above is actually in powers of rj where 77 is our small parameter 
and Gjj and £' are linear operators which are defined fully in paper I. For our 
purposes here, it is sufficient to realize that is actually a first order differ- 
ential operator with respect to the velocity phase space variables. Thus, if the 
quasi-linear diabatic solutions are to be considered a leading approximation to 
the solutions of the master equation in a perturbation expansion in 17 , it is clear 
that the quasi-linear diabatic solutions should not contain velocity space gradi- 
ents which grow large with decreasing 77. But this kind of behavior is exactly 
what we have discovered here. We have found that as r - « the first derivative 
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with respect to a*, of the time average of the quasi-linear diabatic solutions 
approaches a value which, at = 0 » is proportional to 77-2. Thus, at yu, =0, 

77 Gq O (t 7~*) ^ when operating on the quasi-linear diabatic solutions. 

We must conclude that we have no basis for believing that the power series ex- 
pansion converges in the vicinity of = 0; in fact the evidence seems to indicate 
just the opposite. It seems clear that, in the vicinity of yn = 0, some presently 
unlcnown renormalization of the power series expansion must be constructed. 
However, this renormalization cannot be based upon the Fokker- Planck approxi- 
mation since we have shown here that this approximation does not accurately 
represent the leading term in the power series expansion. We point out the 
possibility, which we feel as a result of this work is quite likely, that the re- 
normalized result may not be Markovian nor governed by a velocity space 
diffusion equation. 
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APPENDIX A 


In this appendix we restate Wiener's generalized Tauberian theorem as 
given by theorem 13b., Chapter V of Widder, but in our notation. We then apply 
Wiener's theorem to Laplace transforms, much in the same way as Widder does 
to obtain his theorem 14 in Chapter V, but with a slight modification which 
allows us to then conclude equation ni.50. 
a). Statement of Theorem 13b 
Let 

!)• gi(^) > 0 gi(^) e L (0 < T < <») 

(functions in L can be represented by absolutely convergent Fourier 
integrals) 


2 ). 


I dr g,(r) 0 (-00 5 X < ®e) 

J 0 


3) . F(t) 2 - c (c finite and positive) 

4) . g 2 ('^) be continuous almost everywhere on (0, «>) 


5). 


6 ). 


E e" < T < ^ ~ 

1=-C« ~ 

bi(p) = P f dr gi(pT) F(t) 

Jo 


exist and be bounded for (0 < p < »). 
Then if h^ (0) exists. 


1 i m 

p -♦ 0 


P f dT gj(p T) F(r) = hj(0) — 

Jo r 


I dr g2(T) 
Jo 


! dr gj(r) 


A.l 
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b). Application to Laplace Transforms 

We introduce the notation, f' (/^, r) = (^, r)/9^ and use the theorem 

stated above to prove the following: If f' (//, r) > - c for 0 1 r < « and is such 
that the integral, 


r dr e-P^ f '(/U, t) 

Jo 


A.2 


converges for all p > 0 (i.e., f’ (/j., t) has a Laplace transform for p > 0), and 
if, 

A.3 


j 

-'o 


dr e-P'^ f (/i., t) 


Pi P 


then. 


- r ds f'(/x, s) ~ A(it) 

^ Jo tt 


A.4 


Proof: 

First notice that. 


^ [ds f'(M, s) = 0 

^ ao ' Jq 


A.5 


and since, 


p f dr e~ (' (fj., t) 

Jo 


1 P 


dr t f ' O, 

Jo 


•)l 


A.6 


then. 


p f dr e-P" = 0 

5 0 Jo 


A.7 
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Thus, if we introduce, 


F(/^.t) = 


(0 < T < 1) 


f'(M.r) (1 < r < ~) 


A. 8 


and complete this proof for F (ai, t), then we will have completed this proof for 
f ' (m , t) as well. 

Following Widder, we chose 


gi('^) - e’ 


and. 


(0 < T < oo) 


(0 < T < 1) 


A.9 


E2(^) ■ 


A. 10 


0 (1 < T < o«) 


Then conditions 1)., 3). and 4). of section a), are obviously satisfied. Condition 
2). holds since 


/,oc 

dr e"'^ = f(l+ix) ^ 0 (-« < x < «) 

0 


A.ll 


Condition 5). becomes 


E .n < 


u. b. 


-1 


EaC^)! = “ 


e - 1 


A.12 


Finally 6). is satisfied since 


hi(p) = P J^dT e-P"^ F(^, r) 
is clearly continuous in the interval (0 < p < ®) and approaches finite limits 


A. 13 
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when p approaches zero or becomes infinite. By equation A.3, hj(p) approaches 
A(yti) as p - 0. For t > 1, and 1 < p < », 


pr .= r + (1 - pr > r + p 


-a 


A. 14 


Therefore, 


Ihi(p)l 1 P [ dr e-P’- |F0 i, r)| 
JQ 


< p e~^P“^^ J* dr e”'^ |F (jj., t)| 


A.15 


~ 0 
Pt 


All conditions of section a), are satisfied and when we apply equation A.l we 
obtain the desired result, equation A.4. This result has been used to obtain 
equation III.50 in section n. 
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APPENDIX B 


In order to construct the H-theorem in section III., in addition to assuming 
a "simple" power spectrum we have assumed that (/x) < 0 where is de- 

fined by equation in.57. In addition we have assumed that )Uj(/x = ± 1) is bounded. 
Since the correlation fimction is normally assumed to have a finite spatial range 
such that its moments exist, it is natural to assume = ±1) is bounded. We 
have not found a clear characterization of the magnetostatic turbulence which 
would guarantee that 5 0 however. In the following we simply list a num- 

ber of expressions for 3Tij(^) with the hope that at least one of these may prove 
useful in future studies of this issue. 

Following immediately from its definition, it is easy to see that. 


^l(M) 



ds s R(s) cos 


0 



B.l 


and therefore. 


^ "" (i^i) 


B.2 


With the introduction of the Fourier transform of the correlation function. 


R(K) = 


1 

= I 

In J 


'\j2n 


e-iKC R(Q 


B.3 


we also find. 
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1 i m - 1 r 

p - 0 V2^ 


clK R(K) 


1 - 2p = 


1 - 2p' 


[p2 + (Km - 1)^] [p^ ^ (Km + 1)^] 


which can be expressed as a Cauchy principle value 


^i(m) = " 


1 / 1 




More explicitly, 


^iCm) - - 


£ 0 y2 7T \M^ 


x' e 


and through integration by parts, 


Jfii(M) = - 


1 i tn J 1 / 

e -* 0 I \fj,2j 


'"‘(JN"* di) 


_L / M rtraHiira) 

VM / ^ 


All of these expressions yield. 


J«i(m) = - 


1 - m" 

(1 + M^)" 
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for the exponential correlation function used in the analysis contained in 
papers I. and U. 


